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Abstract 

We revisit the construction of maximal initial data on compact manifolds in vacuum 
with positive cosmological constant via the conformal method. We discuss, extend and 
apply recent results of Hebey et al. m and Premoselli m which yield existence, non¬ 
existence, (non-)uniqueness and (linearisation-) stability of solutions of the Lichnerowicz 
equation, depending on its coefficients. We then focus on so-called (i, (^-symmetric data 
as “seed manifolds”, and in particular on Bowen-York data on the round hypertorus 
§ 2 x 8 (a slice of Nariai) and on Kerr-de Sitter. In the former case, we clarify the 
bifurcation structure of the axially symmetric solutions of the Lichnerowicz equation 
in terms of the angular momentum as bifurcation parameter, using a combination of 
analytical and numerical techniques. As to the latter example, we show how dynamical 
data can be constructed in a natural way via conformal rescalings of Kerr-de Sitter 
data. 


1 Introduction 

We start with two definitions. 

Definition 1. As Initial Data (ID) (Ad ,gij,Kij ) (i,j,—1,2,3) for vacuum with positive 
cosmological constant A we take a compact 3-dim. Riemannian manifold A 1 with smooth 
metric g\j and smooth second fundamental form K,j which is maximal = 0 and satisfies 

the constraints 

R = KijK ij + 2A ViK ij = 0 (1) 

Here V and R are the covariant derivative and the scalar curvature offjij. 

Definition 2. A Seed Manifold (SM) (Ad, g %31 K tJ ) consists of a compact 3-dim. manifold 
(Ad, g t j) with smooth metric in the positive Yamabe class, and of a smooth trace-free and 
divergence-free tensor Ky on Ad. 

The conformal method is the art of turning a SM into an ID via conformal rescaling [24] , 
In the present setting it remains to be shown that the Lichnerowicz equation 

/ 1 \ I Q 2 

tW:=-(A--fl)0- i A#-^=O (2) 
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has a smooth, strictly positive solution </>, where A and R are the Laplacian and the scalar 
curvature of g VJ , and f2 2 = K 1J K,j. In this case the “physical” quantities 

9ij = <P 4 9i> K iJ = (3) 

indeed satisfy the constraints (fTj). 

In view of the observed small positive value of A, and due to the naturality of the assump¬ 
tion of maximality of the data, we are dealing here with a physically very realistic case of 
the Lichnerowicz equation. It is precisely this case, however, which involves rather intricate 
mathematical problems. Firstly, solutions definitely do not exist for large fl 2 which is rather 
easy to see in principle either from the maximum principle, or by integrating fl2]). On the 
other hand, existence proofs for small fl 2 are subtle, in particular when f2 2 is allowed to have 
zeros [m mi i2Qi mi]- However, in physically meaningful situations fl 2 does have zeros - in 
the axially symmetric (AS) case, on which we focus in this paper and which is simple in other 
respects, O 2 in fact typically vanishes on the axis (cf. Sect. 3). 

There are now available two types of general existence and non-existence results which 
cover the case of present interest. The first one, due to Hebey et al. [19] (see also mm) 
guarantees existence of solutions if f M Q 2 is small, and proves non-existence if f M Q 5//6 is 
large. In either case, the bounds can be given explicitly in terms of the Yamabe constant of 
A4 and other integrals over A4. However, there is a f2-“gap” which is not covered by these 
results. In the second theorem, due to Premoselli [ST], H is written as H = b£2 0 f° r some 
(fixed) function f2 0 and (variable) constant b > 0, and the result is “gap free”: It is asserted 
that there is a constant 5* G (0, oo) such that ([2]) has at least two positive solutions for all 
fe < &*, a unique solution for b = fe* and no solution for b > b*. Moreover, for every b < 6* 
there is a unique stable, “minimal” solution. We remark, however, that in this theorem there 
is no direct information about 6* in terms of more familiar geometric quantities of Ai. 

In this work we start (in Sect. 2.1) with defining (in Def. 3) (linearisation-) stability 
of solutions of ([2]) and of ID (under conformal deformations), which will be key in what 
follows. In particular we prove Proposition 1 which guarantees instability if Q 2 = K r] K IJ < A. 
Another important issue in our work is “symmetry-inheriting” versus ’’symmetry-breaking” 
of solutions, by which we mean solutions of (|2]) which share (or don’t share) all symmetries 
of the equation. In Sect. 2.2. we prove a simple result (Proposition 2) which ensures 
symmetry inheritance for stable solutions in the case of continuous symmetries. We proceed 
in Sect. 2.3. by reviewing the theorems of Hebey et al. and Premoselli mentioned above. As 
small complements to the latter result, we clarify (in Proposition 3) how the stable, minimal 
solutions of ([ 2 ]) approach zero as b —> 0. Moreover, in Proposition 4 we employ an argument 
from bifurcation theory to show that near the maximal value 5*, there are precisely two 
solutions. 

The core of our paper is Sect. 3 where we apply the results sketched above to certain 
“(f, ^-symmetric data” as introduced and discussed in [T5[ ITS]. There one sets out from an 
AS “twist potential” oo from which there is constructed an AS, symmetric, trace-free and 
divergence free tensor . The SM constructed in this way “rotate” in general, and the 
(Komar-) angular momentum J of any selected 2-surface is given directly in terms of the 
values of oo on the axis via 8ttJ = [a;(0) — (u(7r)]. We focus on two different classes of seed 
manifolds (Ai,gij) as examples: In Sect. 3.3. we consider a “round hypertorus”, i.e. § 2 x S 1 
with a round S 2 . We first review the case without angular momentum where the solutions 
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of ([2]) yield the time symmetric Kottler (Schwarzschild-de Sitter) data. Then we consider a 
Kij of “Bowen-York form” mi as the simplest non-trivial rotating model. Applying the 
results of Hebey et al. [T5] we find (in Theorem 3) that small angular momenta (compared 
to A -1 , and taken w.r.t. the § 2 surfaces) guarantee existence of solutions of (J2J) while large 
ones exclude existence. Finally, we apply Premoselli’s theorem [3Tj. Combined with auxiliary 
results from bifurcation theory, with results on stability and symmetry collected in Sect 2 ., 
as well as with numerical methods, we are able to clarify the bifurcation structure of the 
axially symmetric solutions in terms of the bifurcation parameter b = 3JA/2: Firstly, there is 
a pair of “principal” branches consisting of stable and unstable solutions all of which inherit 
the 0(2) x 0(2)- symmetry of the SM. These branches emanate at J = 0 from the solutions 
0 = 0 and 0 = 1 , respectively, and meet at some marginally stable solution 0 * corresponding 
to a maximal angular momentum J*. Moreover, off certain points on the unstable principal 
branch there bifurcate branches which break the 0 ( 2 )-symmetry along the § direction, and 
which terminate at the Kottler solutions in the limit of vanishing angular momentum. We 
summarize these facts as Conjecture 1, which also includes the hypothesis that there are no 
solutions which break the axial symmetry on § 2 . 

In Sect. 3.4. we consider as SM the standard maximal slice of Kerr-de Sitter. For 
any fixed A, we take a family of (t , (^-symmetric data generated by uj(r, 9, A, J, a, m) = 
Jux(r, 9, A, a, m )/ Jk where 9 and r are “Boyer-Lindquist” coordinates, J is the “true” angu¬ 
lar momentum and ujk is the twist potential generating Kerr-de Sitter with angular momen¬ 
tum Jk = ma( 1 + Aa 2 /3 )~ 2 in terms of its standard parameters m and a. This example is 
particularly well suited to illustrate Premosclli’s result [3TJ: Choosing b = 3JA/2 as above, it 
trivially implies existence of solutions to ([T]) for all J < Jk- More interestingly, it also shows 
that Kerr-de Sitter can be “overspun” in the sense that there exist data with J > Jk as 
long as they remain strictly unstable. The latter is guaranteed in particular by the criterium 
D 2 < A of Sect. 2 . 1 . mentioned above, but in any case for sufficiently small Jk and small 
J — Jk (again compared to A -1 ). These facts are collected in Theorem 4. 

2 Stability, symmetry, existence and non-existence 

2.1 Stability 

In the following discussion we refer to SMs and IDs as defined in Defs. 1. and 2.. As a rule 
the metric and the second fundamental form of IDs will carry tildes. We note, however, that 
the SM in our examples (Sects 3.3. and 3.4.) trivially satisfy the constraints as well, so they 
are IDs on their own. Hence the task here is actually to generate non-trivial IDs from trivial 
ones. 

We first define and discuss here (linearisation-) stability of solutions of ([2]) under conformal 
deformations, which will be crucial in the following results. The linearized operator L$ 
corresponding to ([ 2 ]) applied to some function 7 reads 



3 



Definition 3. 


1. A solution (ft of ^ on a SM (Ml, gij, K, t j) is called strictly stable, stable, marginally 
stable, unstable or strictly unstable if the lowest eigenvalue g in at (ft satisfies g > 0, 
g > 0, g = 0, g < 0, or g < 0, respectively. 

2. ID given by Equ. (J\) are said to have lowest eigenvalue g (under conformal deforma¬ 
tions) if ^ has this lowest eigenvalue at (ft = 1. The ID are called strictly stable if 
(ft = 1 is strictly stable, and analogous definitions for ID apply with the other stability 
properties. 


The natural question raised by these definitions is resolved as follows. 


Lemma 1. A strictly stable solution (ft of (HP on a SM (Ml, g t j, K l3 ) defines via §3) strictly 
stable ID (Ml,(jij, Kij). The same applies to the other stability properties of Definition^ 


Proof. We show, more generally, that only the conformal class of the SM matters for stability 
of the solution of ([2]) and for the resulting ID. We first note that ([2]) is obviously conformally 
invariant in the sense that when <ft solves (J2J) and defines ID (Ml,fiij = (ft 4 gij,Kij = (ft~ 2 Kij), 
then (ft = d~ l (ft solves (J2D on the SM (Ml, fi l j = d 4 , K tJ = r d~ 2 K i fi) and defines the same 
ID. A conformal covariance property also holds for the linearisation operator ([4|) in the sense 
that the rescaling 7 = d _1 7 gives L 77 = d _ 5 L 07 in terms of L^ on Ad. A subtlety now 
arises since the eigenvalue equation L^p = A/i is obviously not conformally invariant (when 
the eigenfunction p is scaled as above) and the same applies to the eigenvalues themselves. 
However, what matters for stability is only the sign (or the vanishing) of the lowest eigenvalue 
<7 To show that this is actually invariant we recall the Rayleigh-Ritz characterisation, 


? = 


inf 

7GC' oo ,7^0 


Im l L <t>l dv 

Im ^ 2 


(5) 


and note that its numerator is invariant, while the denominator is manifestly positive. The 
statement of the Lemma is now obtained by setting d = (ft in the above arguments. □ 


Remarks. 


1. Recalling that the eigenvalues A depend on the conformal scaling of the metric in general, 
we denote by A the eigenvalues w.r.t to the generated ID (i.e. when (ft = 1 in (jl])). 

2. The above definitions of stability under conformal deformations have nothing to do with 
dynamical stability of the solutions evolving from the data. We will return to this issue 
in connection with the Kerr-de Sitter example in Sect.3.4. 

Proposition 1. Let (Ml, fpj, Kifi be ID with volume V and lowest eigenvalue <7 Then 


1 . 



D 2 dT> (A + ?)R. 


2. If LI 2 < A and D 2 ^ A on Ml, then the ID are strictly unstable. 


( 6 ) 
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3. If Q = 0, then q = —A, w/w/e D 2 = A implies q — 0. 

Proof. Combining ([T|) with (J4j) for <f> = 1 we find 

ac = (^ 2 -a-?x (7) 

where ( is the eigenfunction corresponding to the lowest eigenvalue q. As ( has no zeros 
it can be chosen to be positive. Dividing by ( and integrating, we obtain 

[jti 2 -(- 2 \V(\ 2 )dv = (A + q)V ( 8 ) 

Jm 

which implies assertions 1. and 2.. Point 3. is a consequence of the maximum principle 

applied to (0). □ 

Remark. We wish to clarify here an important point which arises by combining Lemma 
|T] with Proposition |T} These results do not imply that any solution of (j3j) with small 
enough on a SM AA leads to unstable ID AT Rather, the tilde on Q 2 in the requirement 
D 2 = f~ 12 Q 2 < A of Proposition [1] must not be overlooked. In fact, the stable examples with 
small but large D (and hence small <f>) will play a key role below. 

2.2 Stability and symmetry 

We recall from the introduction that “symmetry-inheriting” and “symmetry breaking” are the 
properties of solutions of ((2]) of (not) sharing all symmetries of the equation. This behaviour 
is related to (in-)stability of solutions; we will observe it in the Bowen-York example of Sect. 
3.3.3. We note here the following 

Proposition 2. Assume that a SM (A i,gij) and D have a continuous symmetry f, i.e. 

= 0 Pfl = 0. (9) 

where C^ is the Lie derivative. Then all stable solutions of (0|) are invariant as well, i.e. 


= 0. 


( 10 ) 


Proof. We first note that, for any solution 6 of (f2l) . C^<f> is a solution of the linearized equation: 
Using that the Lie derivative commutes with A and R for invariant metrics g tJ gives the 
second equation in 


0 = Ct 




+ + 


D 2 


= L(C^) 


( 11 ) 


while the first one is obvious from the fact that <f solves ()2|) . 

Next, as <f is stable, the lowest eigenvalue of L is non-negative. As already noted in Sect. 
2 . 1 ., the lowest eigenvalue is always non-degenerate and the corresponding eigenfunction has 
no zeros. It follows from (1111) that either = 0 which we want to prove, or that is a 
ground state eigenfunction with eigenvalue zero and without zeros; by changing sign of f> if 
necessary, we thus have 

Li-f > 0 . ( 12 ) 


5 




(13) 

□ 


This we rule out as follows. Since £* is a Killing vector, m can be rewritten as 

) > 0 . 

But as M. is compact, the l.h.s integrates to zero and gives a contradiction. 

We remark that arguments along the lines above can be and have been applied to a large 
class of semilinear and quasilinear elliptic equations on compact manifolds (cf. e.g. Sect.8 of 

0 ). 


2.3 Existence, non-existence and stability 


We adapt here the results of Hebey et al. [18, [191 EQj and Premoselli [3Tj to the present 
context in order to obtain bounds on 9 and its integrals in terms of geometric quantitites, 
which guarantee existence or non-existence of solutions of ([2]). Premoselli’s result also has an 
impact on the relation between stability and symmetry, which we state in Corollary 1. after 
the Theorem. 

We remark that the results [HI . 3T] refer to a more general equation than (|2]) in which R 
and A can be replaced by a large class of functions. A feature of the present equation ([2]) 
already noted in Sect. 2.1. is its conformal invariance which is obvious from the purpose 
which is serves, and which simplifies the (non-)existence criteria. 

Turning now to the results of [T9], we first note that the existence result Thin 3.1 is 
indeed formulated in an invariant way under conformal rescalings of the SM 

provided the test function <p is assigned the conformal weight (p = On the other hand, 

the non-existence result, Thm.2.1. of [T9] is not conformally invariant. We reproduce these 
results below (Thm. [T]) under the simplifying assumption that the SM has constant curvature. 
Needless to say, this restriction breaks conformal invariance. In the examples discussed below, 
the SM S 2 x S of Sect. 3.3. has constant curvature R — 2A while the Kerr de Sitter data of 
Sect 3.4. have not. 

In the existence criterium (fT5lh there enters the Yamabe constant 


Y 


inf 

7eC°°,7^0 


i M mi\ 2 + Ri 2 )dv 


(Im ^ dv ) 


1/3 


(14) 


We note that Y is conformally invariant. Moreover, by virtue of Yamabe’s theorem [27], there 
always exists a scaling which minimizes Y, and such a minimizer has constant curvature. 
Therefore, the definition (fTTD can be reduced to Y = inf (RV 2 ^ 3 ) where V is the volume of 
A4, and the inhmum is taken over all metrics with R = const, within the conformal class. 

We also remark that, in order to have a chance of satisfying ([!]) it is clear that we have 
to set out from SM (A4, g l3 ) which are in the positive Yamabe class, i.e. Y > 0. 

Theorem 1. We take a SM (Ml, gij, K t] ) as defined in Def. (Hj) but with constant scalar 
curvature R. 


1. Assuming that 


Y 6 


J M ^ dV ~ 256A 2 i? 3 K 3 

Equ. (dj) has a smooth, positive solution. 


(15) 


6 




2. Assume that 


[ 05/6 dv RV (16) 

J M 3 5 / 4 A 5 / 6 1 ; 

then has no smooth, positive solution. 

Proof. For fi = 0 the first part is obvious from Yamabe’s theorem [27] , Otherwise, this part 
is a direct application of Thm. 3.1. of [19], observing that the Sobolev constant Sh of this 
theorem is related to the Yamabe constant Y via (Y/8) 3 = 1 /Sh in the present situation, and 
setting the “test function” cp = 1. (Thereby we arc likely to miss the optimal numerical factor 
on the r.h.s. of (fT5jl ). The second part follows readily from Thm. 2.1. of [19] except for the 
fact that fi > 0 was required there. The extension which allows for zeros in is covered by 
Thm. 3 of [TH]. □ 

We now rewrite Premoselli’s results [31 j . 

Theorem 2. We decompose 11 = 6Ho in (G|) (in a non-unique way) in terms of a constant 
b > 0 and a function Ho- The following statements refer to the solubility of (|I|) on a SM 
depending on the choice of b, when Ho is kept fixed: There exists 0 < 6* < oo such that (Q|) 
has 

1. At least two positive solutions for b < b*, at least one of which is strictly stable. More¬ 
over, one of the strictly stable solutions, called 0(6), is “minimal” in the sense that for 
any positive solution 0 ^ 0(6) we have 0 > 0(6). 

2. A unique, positive solution for 6 = 6* which is marginally stable. 

3. No solution for b > 6* 

Proof. These statements just combine Theorem 1.1., Proposition 3.1. (positivity of solutions) 
and Proposition 6.1. (stability) of [31]. (The statement of strict stability in point 1. is not 
explicit in the formulation of the latter Proposition, but contained in its proof). □ 

The following extension of Proposition [2] is an immediate consequence of this theorem. 

Corollary 1. If {M., gif) and H have a discrete symmetry, then the stable “minimal” solution 
of point 1. in Premoselli’s theorem, as well as the unique solution of point 2. of this theorem 
share this symmetry. 

As typical for non-linear equations, we expect bifurcations to occur among the set of 
solutions of ([2]). While the detailed behaviour of this set will depend on g l3 , A and H, 
Premoselli’s theorem indicates that 6 plays a distinguished role as bifurcation parameter. 
The key values of 6 for understanding the structure of the solutions are 6 = 0, and the 
“critical” values by which we mean those for which the linearised operator (III) has a non¬ 
trivial kernel. The latter is the case in particular at 6 = 6*, but in general (and in particular 
in the Bowen-York example in Sect. 3.3.3.) more such critical values will show up. 

We first discuss 6 = 0. While Premoselli’s theorem does not apply to this case, it is known 
that Equ. ([2]) has at least two solutions: 
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0 > 0: We first mention the special case R — 2A where there is the trivial solution 0 = 1; 
however, there are many more (Kottler-) solutions which we revisit in Sect.3.3.1. In the 
general case, Yamabe’s theorem mentioned above guarantees the existence of at least 
one positive solution. We now observe that all regular solutions are necessarily unstable 
in the sense of Def. [3] this follows from point 3. of Proposition [TJ while point 2. shows 
that instability still holds for small hi = bf~ e flo and therefore small b. However, in 
this context it is important to avoid an instructive catch: Premosclli’s theorem asserts 
the existence of at least one stable solution for all small enough b (and therefore, for 
small enough H = M7 0 )- The key to resolving this issue is the same as in the Remark 
after Proposition Q3 namely a tilde: H = 0 _6 fl = &0 _6 fl o . We are led to the conclusion 
that the conformal factor 0 which generates the stable branch of solutions from any 
SM must go to zero when b —> 0 in order to allow fl to violate the instability condition 
fl 2 < A (or its integral). This leads us to the other solution of (|2j) for b — 0, namely 

0 = 0: While useless as conformal factor, the above arguments indicate that this solution 
is the origin of the unique “minimal” branch of Premosclli’s theorem. Proposition [3] 
confirms and clarifies this. 

In the following result the rescaling 0 = &~ 1(/4 0 will be crucial. In terms of this variable, 
we obtain from ([2]) 

W):=-(a-I«)*Pa#-4 = 0. (17) 

More precisely, (TT7j) is equivalent to ([2]) only for b > 0, but we consider the former equation 
for b > 0. Note that the constant b which controlled the size of the momentum term in ([2]) 
now scales the cosmological constant in (TTTlh 

We also introduce the linearisation at some 0, 

Vt = - (a - V) 7 - jA^ 4 7 + (18) 

Proposition 3. For sufficiently small b> 0, the equation [Ftp on a given SM (Ad, Kif) 
has a unique, positive, strictly stable solution 0(6). 

Proof. For b = 0, it can be shown via the sub- and supersolution method [25] that the 
resulting Lichnerowicz equation (1171) has a unique, positive solution 0 O . Next, strict stability 
follows readily from the linearisation 

io 7 = _( A _i fl ) T+ g 7 . (19 , 

In particular, using that the Yarnabe constant of A4 is positive, L 0 has a trivial kernel. This 
allows application of the implicit function theorem and indeed yields the desired conclusion. 

□ 

Remark. For b > 0 we clearly recover here the beginning of the unique strictly stable 
minimal branch of solutions from point 1. of Premoselli’s theorem. Note that, for b —> 0, 
regularity of 0 indeed entails 0 —)■ 0, as anticipated in the discussion above and in the Remark 





after Proposition 1. Bounds on integral norms of </> in terms of 6 ,A, A and fio can be derived 
via Equ. © but will not be given here. 

We now turn to the critical values of b. As the analysis is slightly simpler in terms of the 
variable 0 compared to 0 , we continue working with (TT71) and its linearisation (IIS]) rather 
than with the equivalent original Lichnerowicz equation (]2|) . 

Here the only simple case is the marginally stable (lowest eigenvalue zero) one, which 
arises in particular at the maximal value 6 = 6 *. In this case a simple bifurcation analysis 
leads to the following behaviour of the solutions: 

Proposition 4. Assume that 0 C is a marginally stable solution of O for some value b c . 
Then there is a solution curve (6(s),0(s)) near ijj c which “turns to the left” (i.e. towards 
smaller values of b) at (6 C , 0 C ) . This entails that there is an e > 0 such that for all 6 G 
( 6 * — e, 6 *) there are precisely two solutions, at least one of which is strictly stable. 

Proof. The requirements of the Crandall-Rabinowitz theorem in the form Thm. 3.2. of [13] 
are as follows: 

1 . The kernel v c of the linearisation L c defined in (IIS]) , and of its adjoint, are one¬ 
dimensional at the critical solution if c . 

2 . The derivative d'&/db\b c of the Lichnerowicz operator (ITT]) is not in the range of the 
linearised operator at b c . 

Now 1. follows from the assumption of marginal stability and the fact that L c is self adjoint 
in the present case. Proving 2. is equivalent to showing that 

L d = ^ ( 20 ) 

has no solutions. Assuming the contrary and using the fact that the ground state eigenfunc¬ 
tion v c can be chosen to be positive, we indeed obtain the contradiction 

0 = f 7 L c v c = [ v c L c 7 = j f if 5 c v c > 0. (21) 

J M J M 4 J M 

From the Theorem we conclude that near the bifurcation point ( 6 C , 0 C ) there is a curve of 
solutions ( 6 (s), 0(s)) = ( 6 C + {3(s),ip c + sv c + r(s)) with /3( 0) = /T(0) = r(0) = r'(0) = 0. To 
compute j3"( 0 ) we differentiate Eq. (fTT|) twice with respect to s and evaluate at s = 0 

L c € ~ (7^oC 9 + 56 C A0 C 3 ) v 2 c - MMv* = 0. (22) 

Multiplying this equation by v c , Equ. L c v c = 0 by fj" and subtracting we obtain 

n ... . 4 f M (fl + 56 c ^) vl 

P"(0) = -2AL2- o If- - crcJ c < o. (23) 

Thus, ( 6 C , "0c) is the turning point at which the curve of solutions turns to the left. □ 

Regarding the behaviour of the solution curve near general critical values 6 C (i.e. with 
negative lowest eigenvalue), it depends largely on the precise form of the equation. We proceed 
with discussing examples. 
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3 (£, (p)-symmetric seed manifolds 

3.1 Angular momentum 

We recall here standard material on axial symmetry (AS) and on the angular momentum of 
compact 2-surfaces of spherical topology. A SM (At, gij, Kij) is AS iff the circle group acts 
effectively on At and its set of fixed points is non-empty. This implies the existence of a 
Killing held rji with fixed points along an axis, such that 

C n gij = C rl K i:j = 0. (24) 

The angular momentum J of a compact 2-surface S in an AS SM is given by 

J=^~ [ Kiji/dS'. (25) 

Since our Definition [2] of a SM contains the requirement that K u ' is divergence-free, all 
homologous 2-surfaces have the same angular momentum. This implies that, in order for J 
to be non-zero, the homology group At) must be non-trivial. We also note that the above 
definition of J is conformally invariant. 

When the SM is AS, so are the stable solutions of (J2]) by Proposition [2] above. The same 
then applies to the ID and, by standard ADM evolution, to the evolving spacetime (A/”, g^) 
(/i, v — 0,1, 2, 3). Any AS spacetime satisfies 

R-iiu A g dL^g^ 0 (26) 

where R^ v the Ricci tensor of and rj M is the spacetime Killing vector. 

The angular momentum and its properties can alternatively be discussed in terms of 
spacetime quantities. In particular, the definition (125]) now reads 

J = IT f V'TdSj (27) 

where V /J ' denotes the covariant derivative of 'g^, and dS^ is the volume element of S. 
From (l26]h the integrand of (1321) satisfies 

v„ (v'Y) = -KiT = -Aif. (28) 

We now recover the spacetime version of the invariance result for J: By Gauss’ theorem, 
(1251) implies that all 2-surfaces S which are homologous and bound an AS 3-surface have the 
same angular momentum. 

We next introduce the twist vector oj^ of the Killing vector rj^ 

( 29 ) 

where is totally antisymmetric and ? 0 i 23 = \J \ det g^ \. cu M is curl-free by virtue of (126|) , 
i.e. V[ /t cA] = 0. Hence there exists locally a twist potential cu, defined up to a constant, such 
that uj p, = 
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For an AS 2-surface S C Mi of spherical topology, the twist potential allows the following 
reformulation of the angular momentum (equivalent to (|25]1 and (127|P . 

J = g MAO - w(S)] (30) 

where N and S are the poles of S. 

3.2 (t, (/?)-symmetric seed manifolds 

Bardeen [3j investigated data for rotating stars which, in terms of particle physics terminology, 
enjoy a PT-invariance, i.e. their evolution is invariant under the simultaneous change of time 
and spin direction. Following Dain [15] and Dain et al. [16] who systematically investigated 
such seed manifolds we call them (t, (^-symmetric. This construction can be summarized as 
follows. 

Definition 4. An AS SM (AT g %3 , K i3 ) is called -symmetric (TPSM) if 

1. The axial Killing field rji is hypersurface orthogonal, i.e. e i3 k= 0 where e^i is 
totally antisymmetric and 6123 = \fdetgij. 

2. K lj satisfies 

Kijffrf = 0 and K i; j(f k q jl = 0 (31) 

where q l3 = gij - rffiypj. 

We now state a well-known result which yields an alternative formulation of a TPSM. 

Proposition 5. 

1. Let (At, g t j, K l3 ) be a TPSM. Then there exists a smooth scalar function uj such that 

(a) The axial Killing field rji leaves uj invariant, i.e. rj l X/ico = 0, and 

(b) The extrinsic curvature 

K ij = 4 (32) 

Tj z 

with r) = rj k qk is smooth everywhere, in particular on the axis. 

2. Conversely, let (A4, g t j) be a manifold of positive Yamabe type such that 

(a) (Mi, gij) is AS with hyper surf ace-orthogonal Killing vector rj l . 

(b) There is a smooth function uj which satisfies l.(a), and K^ defined by ffM) satisfies 
l.(b). 

Then (AT g t j, K t j) is a TPSM. 

Proof. Simple calculations and application of the Poincare Lemma, cf. mm ■ ° 
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Remarks. 


1. In contrast to r/* which is hypersurface orthogonal (w.r.t. a foliation of 2-surfaces) by 
definition of TPSM, the spacetime Killing vector rj a which arises from the corresponding 
data is no longer hypersurface orthogonal (w.r.t. a foliation of 3-surfaces) in general. 

2. The twist potential oj was defined in Sect. 13. II for all AS spacetimes. If such spacetimes 
arise from TPSM generated by the scalar function oj via ([32]) . it can be shown that the 
restriction of oj to the initial surface A4 coincides with oj, provided of course that the 
respective additive constants are adapted. This justifies the synonymous notation. 

3. In a coordinate system (p, 9, (p) where p* = d/dip and the metric is diagonal, TPSM 
have K pip and Kq v as only non-vanishing components of K %i . 

For a TPSM we easily obtain from (l32j) 


Q 2 = K i:j K ij 


|Vof 

2 if 


(33) 


This is the key input for the Lichnerowicz equation in the subsequent applications. 

From now on we specify the SM (A4, p^) to have topology 8 2 x §. For all AS SM and ID 
of this topology, we recall from Sect. (13.11) that the angular momentum does not depend on 
the selected § 2 -surface; one can therefore use the term “angular momentum of the SM (ID)” 
instead of § 2 . 


3.3 The round hypertorus. 

Here we restrict ourselves to the metric 

ds 2 = — [da 2 + dO 2 + sin 2 6dip 2 ) (34) 

where a G [0,T] “goes around” the S-direction. This metric obviously has 0(2 ) x 0(3) as 
isometry group. Equ. (I34|) is also the induced metric on a time symmetric slice of the Nariai 
spacetime [29]. We now consider possible choices for oj and on this background. 

3.3.1 D = 0 

We first recall from Sect. 2.3. that, on arbitrary backgrounds (A4, g tJ ), the Lichnerowicz 
equation ([2]) reduces for D = 0 to the Yarnabe equation. This is the equation which minimizes 
the Yarnabe functional (TT4[) . and the corresponding solutions determine data with constant 
scalar curvature R = 2A. For the present background (l34j) with R = 2A, the Yarnabe problem 
has been studied thoroughly [9[ 21, [32j [33]) due to its simplicity, but also due to interesting 
degeneracy properties. We review the results here. 

While (p = 0o = 1 is of course a solution of ([2]), there are in addition k positive solutions 
4>j iff T G (27rfc,27r(fc + 1)] ( k G N 0 , 0 < j < k). For j > 1 these solutions are periodic 
in a with periods P(m, A) = T/j (cf. (137|) below). The resulting metrics (p^gik are known 
as time-symmetric data for the Kottler (Schwarzschild-de Sitter) spacetime which contain j 
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pairs of horizons (each pair consisting of a “cosmological” and a “black hole” one). Explicitly, 
(f) can be determined from the standard form of the 1-parameter family of time-symmetric 
Kottler data via 

~ dv ^ 

ds 2 K = r 2 (-b dd 2 + sin 2 9dif 2 ) = r(a, m , A ) 2 {da 2 + dO 2 + sin 2 9). (35) 

b- 4 /A 

Here m e M, m < 1 /(3x/3), cr = (r 2 — 2mr — Ar 4 /3), and the horizons r b and r c are 
located at the positive zeros of a. Obviously, the coordinates are related via 

— = a -1//2 r(a = 0) = r b . (36) 

dr 

Integrating over the circle and requiring that the period P{m , A) fits on the hypertorus gives 

P(m, A) := 2 / a- 1/2 dr = T/j j 6 N. (37) 

Jr b 

This way the parameter m acquires a dependence on T/j, and therefore the same applies 
to (j) — (j)j. 

All such data have Ricci scalar R = 2A but the manifolds (A i,<p^gi k ) have different 

volumes, and the Yamabe constant Y = inf (RV 2 ^ 3 ) = 2Amin jV^ 3 is “realized” by the 
manifold of minimal volume. For T < 2n, this is necessarily (134|) as 0 = 0o = 1 is the only 
solution of ()2|) . while for all T > 2n the metric <p^gik with minimal volume always turns out 
to be (j)\g ik . 

Note that the Lichnerowicz equation is independent of a while all solutions except for 
4> = 1 do depend on it. In other words, we have here a simple example of “symmetry 
breaking”. In terms of the stability classification Defj3j we find that both (|34|) as well as all 
Kottler data are unstable, in consistency with Propositions [D and [2l As already mentioned 
in Sect. 2.1. the stability classification should be considered as a mathematical tool rather 
than interpreted physically. 


3.3.2 = const. 


An exhausitve analysis of the solutions of the Lichnerowicz equation in this case has recently 
been obtained by Chrusciel and Gicquaud [10]. In particular, it has been pointed out in Thm 
3.1. there that the results of [26] imply that all solutions of ([2]) are 0(3)-symmetric, i.e. they 
only depend on a. 

The assumption = const, leads to interesting problems regarding the bifurcation struc¬ 
ture of solutions. However, it is incompatible with the ( t , (/^-symmetric scheme described in 
Definition [I] on which we focus in this work. To see this, we integrate (fo3l) with the present 
assumptions which gives to = 4J(sin20 — 29) /it, f2 2 = 128T 2 A 3 /7 t 2 and via (fAZi) a second 
fundamental form whose only non-vanishing component is 


= 


8JA 1/2 sind 


7T 


(38) 


This tensor is singular on the axis, however. We therefore consider different choices of ui. 
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3.3.3 Bowen-York data 

The standard setting for Bowen-York data is a flat SM (A 4,gij,Kij) with 

6 

K%j = -gC ki{iJ k n l nj ) (39) 

where n l is a radial unit vector, r is the radius on R 3 , and the constant angular momen¬ 
tum vector Ji = J8 Z i points in the Cartesian z direction HE]. Alternatively this second 
fundamental form can be constructed via (HUD from the AS function 

ui = — 2J(cos 3 9 — 3 cos6>) (40) 


where J is in fact the standard angular momentum, as follows from HUD . We now carry these 
definitions over to the SM (134D on § 2 x S, replacing (HUD by 

Kij = 6A %€ki(iJ k n l nj). (41) 


Here the radial unit vector n l is orthogonal to § 2 (pointing in the “a-direction” in the coor¬ 
dinates ([33D), while the vector J* now reads A 1 / 2 J(cos0sin0, 0). In terms of the construction 
described in Proposition 0, the generating function still reads as above, namely (HUD . 

Inserting in (HUI) shows that the only non-vanishing component of the second fundamental 
form is 

K aip = 3JA 1/2 sin 2 9. (42) 

This tensor is indeed smooth on the axis, for the same reasons which yield smoothness of the 
metric (I34D . 

It now follows from (HUD that 


Q 2 = KijK ij = = 18 J 2 A 3 sin 2 9 = 8b 2 A sin 2 9 

2 T] 2 

where we have set b = 3JA/2. 

Inserting (134D and (HUD into (J2J) we are left with 


(43) 


•m 

A 




(44) 


where 2 A is the Laplacian on the round § 2 . The corresponding linearised operator around 
some (j) reads 


A 


■7 = 


A + B + -0 4 
4 


7b 2 sin 2 9 


7- 


(45) 


We first adapt the (non-)existence result Theorem [I] to this example. We obtain 

Theorem 3. 

1. IfT < 2tt and A\J\ < 0.05, Equ. 1 ) has a smooth, positive solution. The same applies 
ifT > 2n and A|J| < 2.01/T 2 . 
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2. If A\J\ > 0.165, Equ. f44\ ) has no smooth, positive solution. 

Proof. From Theorem Q] via a lengthy calculation, which uses Y = 87t 4 / 3 (T/27t) 2 / 3 for T < 2tt 
and Y > 87T 4 / 3 for T > 2n in point 1. See [3U] for details. □ 

We now apply Premoselli’s theorem and our results on symmetry and stability to Equs. 
dm) and dm - To obtain a complete picture of the bifurcation structure of the solutions we 
have to resort partially to numerical methods. We first determine the “principal” branches 
of solutions which only depend on 9 and which are equatorially symmetric, i.e. we assume 
0 = 0(0) and 0(6*) = 0(7 t — 9). Due to the symmetries of the SM, and by virtue of Proposition 
[21 we know that this class will include the unique stable, minimal branch whose existence 
is guaranteed by Thm. [21 In order to regularise this branch near b = 0 we now adopt the 
substitution 0 = b~ l ^<f introduced already in (fTTlh This gives 


n ; 1 i5 b 2 sin 2 6 n 

B ‘f+Y + ^^ =0 - 


BV>+0> 5 + ^ = o. (46) 


C\J 

-©- 




Figure 1: The stable (green) and the unstable (red) axially symmetric solutions of (l46lh 
(represented by their equatorial values), and the bifurcation points of the secondary branches 
for the sample period T = 5tt (blue). 

As to obtaining the diagram on the left one first shows that, near a pole (9 = 0) there 
is a 1-parameter family of analytic solutions of the form <f{9) = d + d{l — d 4 )9 2 /16 + 0(9 4 ). 
One then “shoots” (numerically) such a local solution towards the equator and adjusts the 
parameter d such that d<f/d9 vanishes at the equator (7r/2). This gives a stable branch 
emanating from 0 = 0 (green in the online version), and an unstable branch (red) emanating 
from 0 = 1 in consistency with Thm. [2] and Proposition [4J Moreover, numerics shows that 
the lowest eigenvalue ? of the linearisation ( 035D with A = 0) decreases monotonically from 0 
to ? = —A (cf. point 3. of Proposition [T]) along the unstable branch. An analogous discussion 
in terms of 0, which makes use of Proposition [3l yields the diagram on the right. 

However, the “principal” branches displayed above cannot comprise all solutions - we 
recall from Sect.3.3.1. that even in the special case b = 0 there is the symmetry-breaking 
Kottler family of data. We therefore now look for solutions of the form 0(a,6*), i.e. we still 
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assume axial symmetry. For small 6 one can in fact infer, via an implicit function-argument, 
the existence of unstable “secondary” branches emanating from each Kottler solution. To see 
where these branches end up, we rewrite an arbitrary eigenvalue A of the linearised operator 
flT3) at some solution (f)(6) on the unstable principal branch in terms of the eigenvalue A g of 
the truncated operator f (l45l) for A = 0). With a corresponding separation of variables in the 
eigenfunction fi = p(a,9 ) = p a (a)pe(9) we find 


IW-iV+k 

A \Tj J A 


j = 0,1,2... 


(47) 


Bifurcation theory (in essence the arguments of Proposition [2]) now shows that symmetry 
breaking can only occur at solutions where the linearised operator has a zero mode. Setting 
now A = 0 and using the numerical observation that A g changes monotonically from 0 at 
6 = 6* to some negative value (depending on T) at 6 = 0 we find that ( 1471 ) has in fact 
k solutions j = 1,2when T £ (27rk,27i(k + 1)] (Fig.l. shows the bifurcation points 
corresponding to j = 1 and j = 2 in the case T = 57t). Further numerical calculations (which 
are interesting on their own and will be described in detail separately |5j) now in fact reveal 
that from the bifurcation point labelled j (1 < j < k ) on the unstable principal branch, there 
emanates a secondary branch of solutions of the form f>j(a, 6) which continues till 6 = 0 and 
terminates at the Kottler solution with precisely j pairs of horizons (cosmological and black 
hole). 

We remark that we cannot rule out the existence of non-axially symmetric solutions. In 
particular, Thm. 3.1. of Pi has no obvious extension to the present case. However, the 
existence of such solutions is unlikely as we do not find corresponding bifurcation points on 
any of the known branches. 

We summarize the above exposition in the following Conjecture. The status of point 3. is 
“truly conjectural” in the sense just mentioned, while points 1. and 2. are “facts”. However, 
we wish to reserve the term “theorem” to a forthcoming publication [5j where we hope to 
present a complete analytic proof of existence of the solutions, but in any case a detailed 
numerical analysis. 

Conjecture 1. The smooth, positive solutions of the Lichnerowicz equation 0 have the 
following properties: 

1. There exists 6* (6* ~ 0.238,) such that for 6 £ (0,6*) there is precisely one stable solu¬ 
tion and one unstable solution which only depend on 6 and are equatorially symmetric 
(the stable and the unstable principal branch). Moreover, for b = 6*, there is a unique 
marginally stable solution with the same symmetries, while for b > 6* there is no solu¬ 
tion. In the limit 6 —)• 0, the stable solutions tend to zero like 6 1//4 . 

2. IfT £ (2nk, 27r(k + 1)], there exist k values bi...b k , with 6* > bi > b 2 > ... > b k > 0 such 
that, from each point on the unstable principal branch corresponding to bj, there bifur¬ 
cates a branch of unstable solutions depending only on a and 9 (secondary branches). 
Each such branch continues till 6 = 0, where its end point represents the Kottler solution 
with j pairs of horizons. 

3. All solutions are axially symmetric (i.e. independent of p>) 
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We finally note that the maximal value fe* ~ 0.238 obtained from numerics (cf. Fig. 1.), 
corresponding to A| J\ ~ 0.159, significantly exceeds the existence limit given in Thm. 3.1., 
while it comes remarkably close to the non-existence limit, Thm. 3.2.. 


3.4 Kerr-de Sitter 

The Kerr-de Sitter (KdS) 4-metric reads, 


a 


ds A = g^dx^x" = —^ ( dt 


a sin 2 9 , \ p 2 2 , P 2 , o2 


K 


dp ] T —dr T — d9 T 


a 


X 


X sin 2 0 
P 2 


adt 


r 2 + a 2 


K 


dp 


( 48) 

in terms of “Boyer-Lindquist”-coordinates with constants m, a, and k — 1 + Aa 2 /3, and in 
terms of the functions 

Ar 2 cos 2 0 

a = (r 2 + a 2 )(l-) — 2 mr, p 2 = r 2 + a 2 cos 2 6, y = 1 H-(49) 


where a generalises the synonymous function in 3.3.1. 

The constants m and a satisfy bounds given by the extreme solutions, but also “absolute’' 
bounds in terms of A alone, cf. [DEI El E] for a discussion. 

The twist scalar of (HHj) as defined after Equ (1291) takes the form 


uk = ( cos 3 6 — 3 cos 9 — 


a 2 cos 9 sin 4 9 
P 2 


(50) 


where 

ma 
Jk = —y 

K z 

is the angular momentum as defined in Sect. 13.11 The bounds on m and a entail a 
bound on Jk, again in terms of the angular momenta Je of the 1-parameter family of the 
extreme solutions; the latter satisfy the absolute bound A J E < A J max ~ 0.170 saturated for 
one particular extreme solution. Interestingly, this value exceeds the non-existence bound of 
Thm. 3.2. (valid for TPSM solutions of ifTTli on the round hypertorus (134D ). 

We now recall how (148|) arises from ( t , (^-symmetric data. The slice t = const is maximal 
with induced metric 



j 2 k j iij 2 (dr 2 . d9 2 \ sin 2 9 
ds K = gp:dx dx J = p -1-+ 


V a 


X J K 2 p 2 


X ( r2 + a 2 ) 2 — era 2 sin 2 9 dp 2 . (52) 


This metric still fits on a 3-manifold A4 of topology S 2 x S D2j which can be seen as 
in the Kottler case. As before we restrict ourselves to the region r G [rb,r c ] where cr > 0, 
which is bounded by the black hole- and the cosmological horizon. At a = 0 the metric (152|) 
is regularised by replacing r by a defined via (136|) but with a from (149|) . Thus the metric 
becomes periodic in a with period 


P( A, m , a) 


a 2 dr. 


'r b 


(53) 
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As we are not interested in the complete set of solutions here, we take P equal to the circum¬ 
ference T of Ad which leaves us with just one pair of horizons. 

From Proposition 0 and using Remark 2. after this Proposition, we can now determine 
the second fundamental form Kf( via co K = uk\m 

With these preparations we now construct new data as follows. 

Definition 5. We set out from non-extreme KdS data (Ml, gjfj, Kfj) as SM with angular 
momentum J K and with Kjf generated via TWt) from oj k . For some J G R we define cu(J) = 
Ju>k/ Jk and Kij(J) again via HM) . 

ft follows that Kij(J) = KfjJ/Jx and Tl 2 (J) = Kij(J)K l i (J) = J 2 Al 2 K /J\. Note that 
(I5TT) does not hold for J Jk- 

Applying Premosclli’s theorem (Thm. [2]) and Propositions 1, 3 and 4 now yields 
Theorem 4. In the setting described in Definition 0 we claim 

1. There exists J* > Jk such that, for J e (0, J*), Equ. (U|) with Q(J) has at least two 
positive solutions, one of which is minimal and stable. Moreover, for J = A there 
is a unique marginally stable solution, for J slightly below J* there are precisely two 
solutions, and for J > J* there is no solution. In the limit J —> 0, the fainily of 
minimal, stable solutions tends to zero like (AJ) 1 / 4 . 

2. If 

[ n 2 K dv < A.V k (54) 

JM 

for KdS data with Jk, FLk and volume Vk, it follows that J* > Jk- 

3. Inequality J5f\ ) (and therefore the conclusion of point 2.) always holds for sufficiently 
small Jk- 

Proof. We first note that Kij(J) is smooth for all J 6 R. Point 1. follows now trivially 
from the results stated before the theorem and makes no direct reference to the properties 
of KdS. We have stated this point explicitly to illustrate how Thm. fallows to deduce the 
existence of a large family of solutions from a single one. Regarding the case J = 0 we 
note that, apart from 0 = 0, non-trivial solutions definitely exist as well, namely solutions 
to the Yarnabe problem [27]; however, we have no information about their multiplicity and 
propagation to J ^ 0 here. To prove 2. we conclude indirectly: Assume that J* = Jk 
which means that, within the 1-parameter family of Afij (J)-tensors generated from oj(J) with 
J 6 1 as in Definition [5] the KdS tensor Kff had in fact the maximal angular momentum 
permitted by Premosclli’s theorem. Then 2. of that theorem would imply that the KdS data 
are marginally stable. However, (1541) together with Proposition [T] implies strict instability, a 
contradiction. (Note that this Proposition applies here directly since KdS are ID rather than 
just a SM). 

For the final point 3. it suffices to show that f l 2 K < A for small Jk- This is intuitively 
clear as Q 2 K is of order J\ near Jk = 0. To see this in detail, it is useful to rescale all variables 
and constants to the dimensionless quantities 

f = \/A r 6 = 6 <p = ip, a = VX a fh = \fAm Jk = A Jk (55) 
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In terms of these variables, the KdS metric and the terms characterising rotation take the 
forms 

ds 2 K = A ~ l ds 2 K luk = A -1 a Jr — A £l 2 K (56) 

where all quantities on the l.h.s. are functions of A ,m,a,r and 9, while all quantities with 
bars can be written in terms of the functions f, 9 and constants fh and a only (and hence do 
not depend explicitly on A). 

Using (1331) we now observe that Qk can be written as 9, fh , a) = Jk&k, o(f, 9 , fh, a). 
We first show that Vt 2 K 0 is analytic in all arguments. This follows from the fact the denomi¬ 
nator in (1331) only vanishes on the axis where, however, it is “regularized” by the zeros of the 
numerator in the same manner as in the Bowen-York example, and we still have Cl 2 K = f sin 2 9 
near the axis with an analytic function /. Next we recall that, for regular KdS data, fh and |a| 
are bounded from above (by a certain number). Therefore, Vl 2 K 0 (f, 9, fh , a), being an analytic 
function on a compact domain, is bounded from above (by a number). This implies that 
VLk — Jk&k ,o can be made as small as needed by decreasing J K . 

We conclude that the KdS metric can indeed be “overspun” in the sense that, for small 
Jk, one can put more angular momentum than (13T]) on the background (1521) . (We remark 
that our notion of “overspinning” has nothing to do with attempts of exceeding the angular 
momentum limit of extreme Kerr black holes, cf. d]). □ 

We finish with two remarks. 

“Conformally relaxed” Kerr-de Sitter data. Premoselli’s theorem and the previous 
one imply that for any KdS seed data with Jr small enough, there exist stable, mini¬ 
mal data which are conformal to these KdS data. We call them “conformally relaxed”, 
(while the unstable KdS data themselves are “conformally excited”). In view of their 
stability, they will necessarily (by Proposition [2]) inherit the axisymmetry of the KdS 
seed, and they will have the same angular momentum by virtue of the conformal invari¬ 
ance of (j25i) . These data will very likely be non-stationary - in any case, again due to 
their stability, they cannot be a member of the KdS family with different parameters, 
as follows from the proof of the above theorem. In case the data are really dynamic, 
it would be interesting to determine their time evolution. Due to axisymmetry, this 
evolution would preserve the angular momentum. Therefore, it is not impossible that 
the resulting spacetime settles down to the same member of the KdS family one started 
with. 

Data with marginally outer trapped surfaces. Within the conformal method, there 
has been considered the problem of finding ID which contain (marginally) outer trapped 
surfaces ((M)OTS). This was accomplished by imposing suitable boundary conditions 
on the SM, first in the asymptotically flat context [HI Hi but recently also for compact 
manifolds with boundary [22J [23J. However, this work does not cover the Lichnerow- 
icz equation with the present sign of the coefficient of (j ) 5 . It would be interesting to 
handle this case as well [133] . This seems to require combining the techniques of the 
aforementioned papers with those of Hebey et al. [19] and/or Premoselli [31] . 

In the examples considered in this paper, both the Kottlcr as well as the KdS family of 
data contain minimal surfaces. In the former case, the minimal surfaces are expected 
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to turn into MOTS when a small angular momentum is added. In the same manner, 
MOTS should arise when the angular momentum of KdS is slightly reduced or increased 
beyond the stationary value as described above. It would be interesting to settle this 
in the generic case. 
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